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We present a systematic study of the leading twist structure functions of the deuteron, 
F-F J and in a fuUy relativistic approach. Our study is based on a reahstic Bethe- 
Salpeter amphtude for the deuteron, which is obtained as a solution to the homogeneous 
Bethe-Salpeter equation with a realistic A''A'^ kernel. Particular effort is made to connect the 
^ ' structure functions to the densities of the appropriate charges and currents. This allows for a 

systematic comparison between the relativistic and nonrelativistic calculations, by analysing 

10 ■ 

. the same densities in both approaches. Thus, the sources of the relativistic effects in the 

00 , structure functions are understood and clearly distinguished from variations caused by the 

\ differences in the model parameters. We present both the formalism and extensive numerical 

0^ ' calculations for all steps of our analysis. 

' We find that the nonrelativistic and relativistic calculations are qualitatively very much 

alike. However, three main features systematically distinguish a consistent relativistic ap- 
^H. proach from the nonrelativistic one: (i) the binding effects are larger, (ii) the effect of Fermi 

, motion at high x is stronger and (iii) the relativistic description of the structure functions 

fef^2 is fully consistent, unlike the nonrelativistic approach, which is internally inconsistent 
and violates the fundamental sum rules. 



X 
J3 



I. INTRODUCTION 

This paper is the second of the two papers devoted to a study of the deep inelastic lepton scattering 
(inclusive electroproduction) on the deuteron within the Bethe-Salpeter (BS) formalism. The first paper ||^, 
presented a formal approach to the deep inelastic scattering within the BS formalism. The emphasize was 
on the self-consistency of the method, and the development of all aspects of the formalism so that it could 
be applied to a study of the realistic cases of reactions. Basic analysis was performed utilizing the operator 
product expansion method in the leading twist approximation. Subsequently, the results were obtained in 
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the space of the moments of the structure functions. Formal consistency of the approach was estabhshed and 
important sum rules were proven. However, all numerical estimates were performed in the "scalar deuteron" 
model which is rather far from reality. At the same time, it was stated that the following tasks have to 
be done in order to complete the study. First, explicit formulas for the deuteron spin- averaged and spin- 
dependent structure functions have to be derived, and, second, the structure functions have to be computed 
in a realistic meson- nucleon theory. The present paper deals with both of the tasks. 

Utilizing the realistic BS amplitude of the deuteron, we calculate all the leading twist structure functions of 
the deuteron, and analyze the manifestations of the relativistic effects in these structure functions. Some of 
the preliminary results, including solution of the BS equation with realistic iVA^-interactions, were discussed 
previously |^-^. However, this paper is not just a simple sum of results reported elsewhere, nor is it 
based solely on the direct application of the formalism developed in ref. but rather presents a better 
understanding of the formal and physical aspects of the process. 

It is imperative to mention now that in recent years a number of studies of the deep inelastic scattering 
of leptons off the deuteron have been done, referring to various manifestations of the relativistic features 
of the deuteron P-pX|]. These studies certainly influenced our analysis during extensive studies in the time 
following the first publication Next, we found it convenient to use the formalism developed in refs. 
in the late 1970-th (see also [Q). This formalism is more convenient if we deal with the structure functions, 
rather than the moments of structure functions (although it is mathematically equivalent, of course). Indeed, 
the formalism of ref. 0| allows us to calculate the physical moments of the structure functions using the 
deuteron amplitude under the Wick rotation. However, the numerical transform from the moments to the 
structure functions (the inverse Mellin transform) is not a well-defined operation. The new formalism allows 
for numerical treatment of the "inverse Wick rotation" of the Bethe-Salpeter amplitude of the deuteron for 
an application to the calculation of the structure functions. 

Another new distinguishing feature of the present paper is the special accent on the systematic comparison 
of the relativistic and non-relativistic approaches. To do this, we consider a connection of the structure 
functions to the densities of appropriate charges and currents, and then analyze these densities in both 
relativistic and nonrelativistic approaches. In this way, we are able to trace the origins of the relativistic 
effects and demonstrate inconsistencies in the nonrelativistic approximations. The paper contains a number 
of illustrations for every step of our analysis. 

Finally, results for the structure functions published previously have been recalculated using the new 
method and some of them have been corrected slightly. 

The paper is organized as follows. In Section ||, we define the densities of the currents which are used later 
in the analysis of the deep inelastic scattering. For this purpose, the realistic BS amplitude for the deuteron 
is presented, and the relativistic and nonrelativistic expressions for densities are defined and calculated 
within the realistic models. In Section 111 the relativistic formalism of the approach to the deep inelastic 
scattering on the deuteron is developed, including explicit relativistic and nonrelativistic formulae for the 
structure functions, and the sum rules for the structure functions are analyzed. Section ^ contains results 
of the calculations of the structure functions , 2 and gf* in the realistic models, both relativistic and 
nonrelativistic. In Section ^ the summary of results of the paper is presented. Two Appendices contain 
important technical details. 
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II. ANATOMY OF THE DEUTERON 



A. Realistic A'^A'^-interaction and Bethe-Salpeter amplitude 

The basis of our approach to the relativistic description of the deep inelastic scattering off the deuteron is 
a nucleon-nucleon Bethe-Salpeter amplitude. General self-consistency of the approach has been analysed in 
our previous paper The way to construct relevant matrix elements has been shown and important sum 
rules have been proved, both without any reference to the particular model for the BS amplitude. A naive 
numerical estimate has been done within the scalar deuteron model which is found to be in a qualitative 
agreement with the nonrelativistic theory of the deep inelastic scattering [p"5l-[l7|] . 

Apparently, for the realistic calculations of the deuteron structure functions we need the realistic BS 
amplitude. The realistic BS amplitude of the deuteron is such an amplitude which provides us with a good 
description of bulk of the deuteron properties. This is the same ideology as it is in constructing the realistic 
wave function of the deuteron |]l8[-^. The most consistent way to obtain the realistic BS amplitude is to 
solve dynamical problem within the realistic model of the A'^A'^-interaction. The realistic A^iV-interaction 
still cannot be derived from the underlying fundamental theory, QCD, or its nonrelativistic approximations. 



such as the chiral perturbation theory. Although some progress has been achieved in recent years 23 24[|. 
Alternatively, the parameters of the A^A^-interaction can be fixed from the available experimental data on 
nucleon-nucleon scattering (phase shifts analysis), if the general form of interaction is somehow fixed. This 



is a traditional way of action [y_8 19 25|J2^. Unfortunately, the parametrization of the A^A^- interaction is 
dependent on the choice of a dynamical equation. In spite of certain similarity of the sets of parameters in 
all realistic models there is no universal potential for all nonrelativistic and relativistic potentials. 

Only one parametrization of the A'^A'^-interaction is available for the BS equation, the parametrization of 
Fleischer and Tjon j25j. This parametrization is probably already due for a revision |2^, incorporating new 
data on the nucleon-nucleon phase shift analysis |p7[| . However, this parametrization has been used as the 
basis for our recent calculation of the BS amplitude of the deuteron |Q. The meson parameters (masses, 
coupling constants, cut-off parameters) have been taken to be the same as in ref. [^5[^, except for the 
coupling constant of the scalar (T-meson, which has been adjusted to provide a numerical solution of the 
homogeneous BS equation. The chosen set of meson parameters requires some minor adjustment, in view 
of two circumstances. First, here we use a simplified form of the propagator of the vector mesons, omitting 
k^k^ / H^-levTB.. Second, the different numerical procedures in solving the eigenvalue problem for the BS 
equation can also affect the value of parameters. 

A detailed formalism preceding the numerical solution of the BS equation for the deuteron is presented 
in ref. The Fredholm system of "Wick rotated" equations with all meson exchanges is solved by an 
iteration procedure with two dimensional gaussian integration (see e.g. ref. [p8|). Eight components of the 
BS amplitude of the deuteron: 

■0pl, V'al, i^vl, IpaO, tpa2, ^tl, V'tO, tpt2- (2.1) 

have been computed as scalar functions of two variables, po and |p|. Now these components are available in 
the form of an analytical parametrization p9[ . 
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Since the one-boson exchange potential |^,^ was used only with a minor adjustment, our solution does 
not contain new physics (or different physics) than the pioneering paper. All parameters are presented 
in Table 1. The set of parameters includes the nucleon mass, m, the deuteron binding energy, ed, and, 
therefore, the deuteron mass, Md = 2m + en- The coupling constants are shown in accordance with our 
definition of the meson- nucleon form- factors: 

= (2.2) 

where A is a cut-off parameter (see Table I). 



B. Static properties and densities 

In order to compare different models of the deuteron and define to what degree each of the models is 
realistic, it is natural to calculate various observables within these models and compare them with each 
other and with the experimental data. The canonical way to do this is to calculate the nucleon contributions 
to static observables, such as mean square radius, magnetic and quadrupole moments, etc. p^- |2^ , |30[ | . (The 
mesonic corrections are quite small and depend upon additional model assumptions.) Most of the modern 
realistic models give quite good description of the static properties, or have a decent explanation in the case 
of a minor discrepancy. The observables of the reactions, such as formfactors and structure functions, are 
more model dependent. Those are often used as a tool for discriminating between models or certain ways of 
using models for calculations. 

A direct comparison of the wave functions and amplitudes usually seems to be less meaningful. Still, in 
many cases particular properties of the wave functions and amplitudes directly indicate what will happen 
when observables are calculated. A good example here is the D-wave admixture, which directly affects 
values of the magnetic and quadrupole moments, tensor analysing power, T2Q, and spin-dependent structure 
functions of the deuteron, and bf'. Another example, the high momentum "tail" of the wave function 
(amplitude) dominates in the observables for some kinematic conditions with a high momentum transfer. 
However, the interpretation and properties of BS amplitudes are very different from the ones of the wave 
functions, a direct comparison between them is impossible. 

Recently, it has been argued that a new intuition can be developed in working with the BS ampli- 
tude 1^,^. The approach is based on dealing with the charge and current densities. The same densities 
also can be calculated in the nonrelativistic approach with wave functions. Thus we gain a possibility to 
compare the models by comparing densities. Besides densities are more directly linked to the observables 
than wave functions and amplitudes. Following this guiding idea, we calculate the charge densities which 
are related to the various structure functions in deep inelastic scattering. 

The nucleon contribution to the deuteron observable, (O), in many important cases is defined by the 
triangle diagram (see Fig. |I]): 

^^Tr {*A,/(po, p)d{pi,q, fc)*M(po + ko,p + k){p2 - m)} , (2.3) 

where pi^2 — Pd/2±p — (A/d/2 ± po, ±p), Pd = {Md,0) is the deuteron momentum in the rest frame, 
p — {po, p) is the relative momentum of nucleons, 0(jJi, q, k) is an appropriate operator and ^'m(poi p) is the 
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Bethe-Salpeter amplitude for the deuteron with M being the deuteron's total angular momentum projection 
(see Ref. for definition and conventions for 5'm(po,p))- The operator 0{pi,q,k) is written in a general 
form, depending on nucleon momentum, pi and two external momenta, q which does not sneak into the 
lower part of the diagram and k which adds to nucleon momentum. In the present paper we consider only 
the case with fc = 0. 

The explicit form of operators and structure of the deuteron matrix elements relevant to the deep inelastic 
scattering, were studied previously |p^ , p7| ,p|-[7|Jl[| . One of the methods to find a form of the operator for deep 
inelastic scattering is to adopt the Wilson's Operator Product Expansion ||3^ and use it within the effective 
meson nucleon theory . Other possible approach is based on a parametrization of the operators 

in the most general form with analysis of all structures phenomenologically or within models for quark- 
nucleon amplitudes |p|-f7|, p^|j33| . Both approaches lead to the same results in "convolution approximation" 
which assumes that the deformation of the nucleon operator off-mass-shell can be neglected. Very small and 
strongly model dependent corrections to the convolution are qualitatively similar in both approaches. In the 
present paper we do not concentrate on the differences between these two approaches, i.e. our analysis is 
not going beyond the convolution approximation. 

Thus, the relevant operators are vector, 7^, and axial vector, 757^^. Let us define the matrix elements of 
the following components of these operators: 



(7o)j\/ 

(73>J\/ 

(7570) J\/ 

^ (7573)j\/ ^ 



> = 



2M 



D 



Tr{^'Af(po,p) 



70 
73 
7570 
L 7573 ) 



>*m(po,p)(P2 ~m)} , 



(2.4) 



where notation in the l.h.s. are obvious and we use this notation in what follows. There are three main 
combinations of these matrix elements which are important for deep inelastic scattering and which are 
discussed in the present paper. They are presented below together with the correspondent structure functions 
of the deuteron: 



9? 



{(7o>A/=i + (73)a/=i} - {(7o)Af=o + (73)Af=o} 

{(7570) A/=l + (7573)Af=l} , 



(2.5) 

(2.6) 
(2.7) 



where — > means "corresponding to" and does not have a direct mathematical interpretation. Explicit con- 



nection of these matrix elements to the structure functions will be discussed in Section III. The important 
circumstance now is that these matrix elements present charges and currents, vector and axial vector, cal- 
culated for the deuteron states with different combinations of the total angular momentum projections, M . 
Properties of these matrix elements, including rigorous sum rules for the charges, are used later to discuss 
the structure functions. 



Let us start from the matrix elements involved in eq. (2.5). The first term, (7o)Af, is the charge of the 
conserved vector current and since it is independent of M the following sum rules can be immediately written 
down (see e.g. p^Jl|): 
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^^(7o)m = (7o)a/ = 1. (2.8) 

M 

Obviously this sum rule is also valid in the nonrelativistic approaches, since it is just a normalization condition 
for both the wave functions and BS amplitudes. To make more meaningful use of this matrix element and 
compare our approach to other realistic approaches, we define charge density in the BS formalism as (see 
Appendix ^ for explicit form): 

(7o>''^(p) = T^-4]n-E / y^Tr{*M(Po,p)7o*M(Po,p)(P2-m)}, (2.9) 

>^^(p) = l, (2.10) 
which in the nonrelativistic limit corresponds to the momentum density [|r; 

(7o)^«(p) = ^ E I^m(p)|' = (^.'(IpI) + ^'(IpD) , (2.11) 

M 



1 

An 



dp{lor'\p) = 1, (2.12) 



where 5'm(p) is a nonrelativistic wave function of the deuteron (do not confuse with the Bethe-Salpeter 
amplitude, ^'m(poi p)0 with u and w being its S- and D-wave components. The densities, (7o)(p), calculated 
in three realistic models |l8| , ^ p| are shown in Fig. || as functions of |p|. Certain model differences at 
IpI > 0.5 GeV are present, however there is no distinguishing feature of the relativistic density. 

Similarly, we consider other components of the vector and axial matrix elements. The relativistic den- 
sities, (73)^'^ (p), (7570)^'^ (p), etc, are defined by replacing 70 with 73, 7570, etc, and taking appropriate 
combination of the matrix elements with different M in eq. (2.9). The explicit expressions for all densities 
in terms of components of the BS amplitude are presented in the Appendix 

The nonrelativistic densities are calculated using methods of nonrelativistic reduction ||l^,|l|,|| • Here we 
present the most interesting of them. A remarkable feature of the non-relativistic densities is that they have 
the same angular dependences as the corresponding relativistic densities. 

The 3rd spatial component of the vector current in the nonrelativistic limit gives: 

(73)"^(P) -/l-i: l*M(P)r + o(^-)^ (u^(lpl) + , (2.13) 

where "d is polar angle of vector p. Note, this density has a suppression factor of ^ |p|/m-, due to a mixing of 
upper and lower components of the Dirac spinors by non-diagonal matrix 73. Since the angular dependence 
of the explicit form of the relativistic density (73) '^'^(p) is also absorbed in the factor cost? (see Appendix 
^), the integrals of both relativistic and nonrelativistic densities (73) over dp are zero: 

1 

j dp(73)^^(p) = j dp(73)'^''(p) = (73) « j d{cosd)cosd = 0, (2.14) 

-1 

These densities for realistic models are show in Fig. ||at 6* = 0. Basically they just reflect the behavior of 
the charge densities from Fig. |^ in accordance with the nonrelativistic formula (2.13). For illustration, the 
curve representing the BS density from Fig. |^ multiplied by the factor |p|/r7i is also shown (dash-dotted). 
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Surprisingly enough, it barely can be distinguished from the exact (73)'^'^(p) (solid line) even at momenta 
higher than m. 

Two examples of spin-dependent densities: 

{lo)Zli{p) - (7o)fA(p) ^ ~lP2{cos^)wi\p\) (2V2u{\p\) + w{\p\)) , (2.15) 

(7573)f/^i(p) ^ u'i\p\) ^w'i\p\) + P,icos^)wi\p\) [w{\p\) - V2u{\p\)) , (2.16) 
where ^2(2;) is the Legcndre polynomial. We also can easily write down sum rules for the spin-dependent 



densities ( |2.15| ) and ( p.iq ): 

rfp{(7o)j^/^i(p) - (7o)j^A(p)} (X / d{cos^)P2{cos^) = 0, (2.17) 



^ / dp(7573)^^(p) = 1 - IwD. (2.18) 



where wd is the weight of the D-wave in the deuteron. The relativistic analogue of the sum rule (2.18) can 
be used for an estimate of the "admixture" of the D-wave in the relativistic formalism which otherwise does 
not allow for probabilistic interpretation. Numerically we have: 

^ J dp (7573) ''^(P) ^ 0.9215, (2.19) 

which gives us an estimate of wd ~ 5%]^ 

The realistic model densities ( |2.15| ) and ( ^.16| ) are presented in Figs. ^ and ||, respectively. These two 
examples confirm the conclusions of the previous illustrations: (i) realistic models are in a reasonable agree- 
ment among each other, providing description of the charge and current densities of the deuteron, and (ii) in 
spite of some model variations at high momentum, there is no distinguishing feature of the densities obtained 
in a relativistic BS formalism. Therefore, we can not expect striking relativistic effects caused by the form 
of the densities. However, this conclusion does not close a possibility for effects generated by the differences 
in the relativistic and non-relativistic description of the deep inelastic reaction. 



III. RELATIVISTIC THEORY OF THE DEEP INELASTIC SCATTERING ON THE DEUTERON 

A. Definitions and kinematics 

We start with the general form of the hadronic tensor of the deuteron with the total angular momentum 
projection, M, keeping only leading twist structure functions: 



^ There can be corrections, Odpp/m^), to the density, {7573)ju=i (p), see refs. ^Jq]. However, their estimated 
contribution to integral (p.lq) is small, < 1%. 



7 



W^Al^ Pd,M) = ( -.9^^ H ^ 1 [XD,Q , M) + - g^^- I Pd,/ - q^ — 



q^ )\ ^ q^ ) Pnq 

^ep.o/3g"5^(M)gf(x,,,Q2), (3.I) 

t^oq 



where q — (y^ 0, 0, — + Q^) is the momentum transfer, — — q^, = ji^Puq) (in the rest frame 
of the deuteron = Q"^ / {iMBv)), Sd{M) is the deuteron spin (see Appendix ^) and F['2 and are 
the deuteron structure functions. Averaged over the M this expression leads to the weh-known form of the 
spin-independent hadronic tensor which is vahd for hadrons with any spin: 

W^iiq, Pd)^\Y1 K'^^I^ Pd,M) (3.2) 



3 M 



It^lA^Df n^^^fp PDq\(„ PDq\ Ff{xD,Q^) 



where F[^2i^D,Q^) are resuh of averaging of the SF F['2{xd, Q^, M). It is easy to find that other structure 
functions can be obtained as a result of other combinations of W^^{q, Pd, M) with different M: 

W;^,{q,PD,M=l)-W;^Jq,PD,M^~l) cx gP{xD,Q^), (3.4) 
W'iiq, Pd,M = 1) - W!^,iq, Pd,M = 0) cx b^JxD,Q^) (3.5) 



The eqs. (3^), (3^) and ( |3.5| ) are the basis for the experimental measurements of the deuteron structure 
functions. However, for the purpose of theoretical studies of the hadronic tensor and structure functions 
the projection technique is more convenient. All relevant formulae for projection technique are presented in 
Appendix More background information on the SF can be found in refs. [p4|-p7[. 

To calculate the hadronic tensor of the deuteron we follow the general formalism of our approach: 

W;^,{q,PD,M) = i J 1*4/(^0, p)W^jI(pi, 9)«'m(po, P)(P2 -m)} , (3.6) 

The nucleon tensor operator, W^^{q,p), has been studied extensively in recent years [[l6y32P j5|,|6|,p|,p^ . We 
use a well established form of the operator, leading to the convolution formula p8[ |: 

W;:L{q,p) = VF{^.}(9,p) + %.i(g,p) (3.7) 

W{,.}iq,p) = ^wl^Aq.p), (3.8) 

W[^.]iq,p) = ^eM-/3g"7''755f (g,p), (3.9) 

where {. . .} and [. . .] denote symmetrization and antisymmetrization of indices, respectively, and gi' {q,p) = 
g^{x,Q^) is the spin-dependent nucleon SF. The hadronic tensor of the nucleon, Wf^^{p,q)^ , is defined as: 

\ ( , 9t^qi'\rpNf , f pq\ ( pq\Fi'{x,Q^) 
w^M^p) = ( -5m'' + ] Pi (a;,Q ) + [Pf, - qf,-^ I \p^ ~ qu-^ I — — , (3.io) 

where x — /{2pq) and F^2 are the nucleon SFs. The small effects of the off-mass-shell deformation of the 
nucleon tensor fippyl are not considered here, since these effects do not affect sum rules for SF and do not 
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noticeably change the absolute values of the SFs. That is why SFs iii ( t^-lC| ) do not depend on p^, but 
and "pq. 

Using projectors (Appendix]^), we extract SFs from the hadronic tensor of the deuteron: 

Q^ M) ^ . / p-F- { Q^) + (3.11) 

J V'T^) VP10+P13 / 2(pio+pi3) 

T?D(^ n-i M\ f '^'^P T^N ( XNTTi ^2^^ Tr{^M(po,p)(7o + 73)^A/(Po,p)(P2-m)} 

i'2 \Xn,Q ,M) = I I ■ ,Q — , (3.12) 

■ [ d^P n( XN-m ^2^ "'"''{^AKPo,P)(7o+73)75^M(po,p)(p2-m)}| 

5i {xN,Q « / TTTTiSi ,Q ^7 ^ 3.13 

J (27r)4 VP10+P13 / 2(pio+pi3) 

where xn ^ Q"^ /I2nii') is the Bjorkcn scaling variable^, i.e. this is x for the on-mass-shell nucleon at rest, 
Pio and P13 are the time and 3-rd components of the struck nucleon momentum. Formulas ( [3.11 ) and (3.12) 



have not been averaged over the projection M since in the present form it gives an understanding of the 



SF hi 2- Indeed, for instance eq. ( 3.12 ) gives two independent "SFs", with M — ±.1 and M = 0, which are 
related to the usual spin- independent SF, F-f , and a new SF, 



F^{XN,Q')^\ E FP{xn,Q\M), (3.14) 

A/=0,±1 

h2{xN.Q^) = F^{x, Q\M = +1) - F^{x, Q^ M = 0), (3.15) 
Fi^{xN, Q^M = +1) = F^{xN, Q^M = -1). (3.16) 

Note, the SF F^{x,Q^,M) is independent of the lepton polarization, therefore, both SFs, F^ and 
can be measured in experiments with an unpolarized lepton beam and polarized deuteron target. In view 



of eq. (3.16), only one of the SFs F^{x,Q'^,M) is needed, in addition to the spin-independent F^{x,Q^), 



in order to obtain b2{x,Q^). The other SF, 6f , is related to the deuteron SF Ff, the same way as fef* is 



related to , viz. via eqs. (3.14), and 62 = 2a;6f 



B. Singularities of the triangle diagram and calculation of structure functions 

It has been shown previously [prj-ITit how a singular structure of the triangle graph (Fig. |l|) rules the 
behavior of the spin- independent SF F^ . In particular, it is found that the relativistic impulse approximation 
satisfies the unitarity and provides the correct kinematical region of the variable xn- However, both these 
properties of the exact covariant amplitude are violated in practical calculations, when nonrelativistic wave 
functions of the deuteron are used. In this case one can refer to the argument that such deviations are small, 
and are not important for phenomenology. At the same time, a realistic Bethe-Salpeter amplitude of the 
deuteron serves ideally for a consistent phenomenological application of the covariant theory of the processes 
on the bound nuclcons. 



^Note that the "native" deuteron variable is xd — {m/MD)xN, however xn is used more often. 
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In order to calculate SFs, ( 3.12 )-( 3.15 ) and analyze the sum rules, the singularities of the triangle diagram 
should be explicitly taken into account. To do that, eq. ( ^.12 ) is rewritten as: 



2Md 
1 



(2^) 



4^2 



N 



PlQ +Pl3 



{Pi — ni^ + ie)'^{p2 — m? + ie) 



-Tr{(/)A/(po,p)(pi +to)(7o + 73)(pi + m)0M(Po, P)(P2 + m)} , (3.17) 



where 0a/(pojP) — {pi ~ '^)'^ m{pot'p){'P2 — fn) is the Bethe-Salpeter vertex function of the dcuteron. 

Analysis of singularities in the complex p2+-plane allows for one analytical integration in (3.17) ||T^. After 
translation into variables which are used in the present paper, this integration is equivalent to picking the 
residue of the second nucleon pole, p20 = ljJ = \/ vn? + or po = — tj, in the complex po-plane when 

both of the following conditions are satisfied: 



< w - p3 < Md- 



(3.18) 



The contribution of the region of p beyond (3.18), to the integral (3.17), is zero, i.e. different poles cancel 
each other. Note, that pio = Md — in the required pole. Calculating residue in ( 3.17 ), one gets: 



1 



2M 



D 



{2n) 



3^2 



xnttl 



Md - lo +pz J 



Q2) e(Afo-cj+p3) 



(3.19) 



rTr{0A/(po,p)(pi +"i)(7o + 73)(pi +"^)0m(po,p)(p2 + m)| 



2wMI{Md - 2w)2 

where the 0-function guaranties the right of conditions ( |3.18 ) is satisfied, and the left condition is always 
satisfied. 

It is useful to rewrite ( 3.19| ) in the convolution form: 



F^{XN,Q^M) = 



(3.20) 



where "the effective distribution" of nucleons in the deuteron is defined by 



fTiy) 



1 



2Md J (27r)3 



5\y 



Md -W+P3 



e(y) 



2ujMfj{MD - 2cj)2 



(3.21) 



Tr {0A/(po,p)(pi + m)(7o + 73)(pi + m)(j>M{po,p){P2 + m)} a/^ 
The SFs i^j^2 ^^^d ^^''^ iio'^ calculated as follows: 

M D /m 



FfixN,Q^) 
b?ixN,Q^) 



dy\ r/^(2/) 
y \ Af^/^iy) 



y 



(3.22) 



FfixN,Q') 



Mn/m 



dy ■ 



/^/^(y) 
Af^/^{y) 



(3.23) 



where distributions f^/^ and A/^/^ are given by 
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Ar/^(2/) = /;'-(2/)-/;'-(y) 



M 
N/D, 



fN/D, 



Similarly, for the SF we get: 



where the effective polarized distribution of nucleons in the deuteron is defined by 



1 



\^\v 



1 



2Md J (27r)3 

Tr {(f>M=i{pQ,p){pi + to) (70 +73)75(^1 + m)(l)M=i{po,p){p2 + m)} ^^^md_^ 



(3.24) 
(3.25) 



(3.26) 



(3.27) 



C. Sum rules for the deuteron structure functions 



Two sum rules can be written for the effective distribution fj^^j 

AIn /rn 



N/D 



1 

47r 



dp{lo)M{p) = (7o)a/ = 1, 







yf^/''iy)dy={D\ieMD)M = l~S_ 



N, 



(3.28) 



(3.29) 



where (Qn)'^ oc iip{x)^fj,d'^ip{x) is the trace of the energy-momentum tensor. Eq. ( 3.2g ) presents the vector 



charge conservation generalized for the deuteron states with different M (see eqs. ( |2.9D -(2.10)). In spite 
of clear physical interpretation, some time ago it was a subject of some controversy [O-H. Indeed, the 



derivation of sum rule (3.28) contains some subtle points and equivalence between it and the expression for 
the charge, (2_^), is a non-trivial fact, particularly, because of the presence of the 0-function in eq. ( 3.21 ). 
This 0-function provides correct kinematics in variable xn but cuts out a part of the integration domain 
in (Pp. This cutting of the integration interval in the polar angle leads to non-zero contribution of the 
matrix element containing 73, which is proportional to cosO. However, the validity of this sum rule has 



been firmly established |ll|-|l3|] . The sum rule ( 3.2£ ) for the first moment of f^/''^ is of a different nature 



it presents the nucleon contribution to the total momentum of the deuteron ||l^,^,|l^,0| and ^at is a part of 
the total momentum carried by the non- nucleon component (mesons). The constant Sn cam not be fixed 
in a model independent fashion, rather it is calculated within any particular model. Self-consistency of the 
theory requires that meson exchange current contribution into the deuteron SF, F.p , exactly compensates 
the loss of energy by nucleons, ( 3.29| ). An important property of sum rules, ( 3.28| ) and ( 3.29|) , is that their 
r.h.s. does not depend on the deuteron spin orientation. 



The sum rules for f^^^{y) and A/^/^(j/) follow from sum rules for f^J^ {y) and definitions (3.24)-(3.25): 
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Mrt /m 



J r/^(y)d2/=i^(7o)M = (7o) = l, 



Mo /m 

I yf'"''{y)dy = \ Y.{d\ {e^r^ \d)m = i-5_ 





Md /m 



Af^/''{y)dy = (7o)a/=i - (7o)m=o = 0, 





Md /m 



yAf''/^iy)dy = {D\ {Q^V \D)m=i {D\ \D)m=o = 0. 



(3.30) 
(3.31) 
(3.32) 
(3.33) 



Sum rules for the deuteron SFs 6f and are the immediate result of combining eqs. (|3|)-(|l3|) and 
(p2|)-(p:23|): 



dxobi (xd) = 0, 



J dxDb^ixo) = 0, 



(3.34) 
(3.35) 



in an agreement with the sum rules suggested by Efremov and Teryaev [ p4[ . 

Sum rule for the spin-dependent distribution relates the integral of the spin-dependent distribution of 
nucleons to the 3-rd component of the axial current, (2.19): 

Mo /m 



/^/^(y)dz/ = (7573)^li. 



(3.36) 



An explicit expression for the distribution function {y) (and therefore of f^^^{y) and /S.f^/^{y)) in 
terms of the components of the Bethe-Salpeter amplitude can be directly established from eqs. ( |3.21 ) and 
(3.27) and formulae for the corresponding densities are given in Appendix H. 



D. Calculation of distributions: Inverse Wick rotation. 



To calculate numerically the effective distribution functions, ( 3.21 ), ( 3.24 ), ( 3.25 ) and ( j3.27| ), we need to 
know the matrix elements over the BS vertex functions as functions of p and at po = Mo/'^ — \Jrn? + 
along real po- Let us remind that the components of BS amplitude have been found along imaginary axis in 
Po-plane (Wick rotation). 

We explain the procedure of calculating the "inverse Wick rotated" matrix element using an example of 



/^/■^(y). Other distributions are calculated the same way. First, let us rewrite eq. ( 3.21 ) in the form: 



r/^(2/) 



d^x 



2Md J (27r) 



Mr 



P3 



e(2/)((7o)fl(p) + (73)fl(p)) 



(3.37) 
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where (7o)pofe(p) ^^"^ (73)pofe(p) ^"^^ densities defined similarly to eq. (|2.9|), but in the "pole approxima- 
tion", i.e. exact integral is replaced by residue at the second nucleon pole. An explicit form of (7o)^ofe(p) 
and (73)pofe(p) is clear from eqs. ( 3.21 ) and ( ^3.37 ). Second, we notice that this density does not have any 



singularities in the complex plane of po. Neither does the matrix element: 

Tr {^m(po,p)(pi +nT-)(7o +73)(pi + m)(j>M{po,p){p2 + m)} 
= {pI - m^ fipl - m2)Tr {*Af(po,p)(7o +73)*Af(po,p)(p2 ~ m)} . (3.38) 



Therefore, to calculate the density in the pole approximation, the matrix element (3.38) can be safely 



expanded into a Taylor series in the variable po around the point po = 0. Third, the coefhcients of the 



Taylor expansion can be calculated, using the known r.h.s. of eq. (3.38) along imaginary po- And last, a 



convergence of the expansion can be checked numerically, by comparing results of calculations up to different 



order in pq. Note, that eq. (3.38) and similar expressions for all other distributions are even functions of 
Po, therefore the Taylor expansion really should be done in p^. In addition, the point pq = should be 
a good point for the Taylor expansion, since in the most physically interesting region of |p|/m <C 1 we 
have |po/"^| ~ \£D/{2m) — |pp/(2m^)| <^ 1. The critical point for any expansion in the nuclear physics is 
usually |p|/m. = 1. However, even at this point jpo/^T-l is still a good parameter for the Taylor expansion, 
\po/m\ ~ |eD/(2m) — (V2 — 1)| ^ 0.4. Parameter |po/"i| is getting close to 1, with |p|/m ^ 1.5. 

The results of calculations of the charge density in the pole approximation, (7o)po;e(p), are presented in 
Fig. ^, where we compare curves for calculations up to ~ Pq , ~ Pq and ~ Pq. As we expected, the procedure 
is nicely convergent up to |p| ^ 1 GeV and with reasonable accuracy can be used up to |p| ~ 1.5 GeV. 
Similar results for the axial density in pole approximation, (7573)po;e(p), are shown in Fig. 0. 

Note, the numerical approximation made in this section, such as limiting number of terms in the Taylor 
expansion, can potentially cause a violation of the exact sum rules. 



E. Nonrelativistic formulae for structure functions 



The nonrelativistic expressions for f^^^iy), Af^/^{y) and f^^^{y) can be obtained by using an analogy 
of the charge densities calculated within the Bethe-Salpeter formalism and corresponding densities calculated 
with wave functions (e.g. (2.11), (2.13), ( ^.15 ) and (2.16)). Actually, the distributions in the BS formalism 



are expressed in terms of densities in the nucleon pole approximation, ( 3.21 ), ( 3.25 ) and ( 3.27 ) and not exact 



densities as (2.9). However, one can hope that the nucleon pole contribution gives the main contribution 
to the density, at least in the non-relativistic region. Such type of assumption is very common for nuclear 
physics (see e.g. 20 j2^,^,^,0). For instance, the well-known result for the spin-independent distribution is 
immediately reproduced (see e.g. [ p"7|Jl^Jl^ ) : 



-(5 y 



(2^)^ 
d^p 

(27r)3 



W +P3 



5[V 



m 

Mp - w +P3 
m 



e(2/){(7o)^«(p) + (73)^^(p)} 
\p\cos6 ^ 



e (y) (1 



-){u'{\p\)+w'{\p\)} 



(3.39) 



The presence of the 0-function on the r.h.s. of eq. (3.39) slightly violates the sum rule (3.3C). However, this 
is not noticeable phenomenologically, since the only region of large momenta, |p| > 0.7 GeV, is affected by 
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the 0-function and it does not contribute much to the norm of the deuteron wave function. We can accept 
such shght effect of the 0-functions, since a nonrelativistic approximation is based on the behef that high 
momenta are not important and if something is wrong there we just ignore it. 
For distribution A/^/-°(j/), we get: 



d\y 



(27r)3 
(H 



P3 



e(y) 



(3.40) 



IpI COS0 



m 



)P2{cos9)^w{\p\){2V2u{\p\)+w{\p\)} 



Again, the sum rule ( |3.32 ) is broken by the presence of the B-function in ( [3.4C ). Neglecting it, one gets 



/ 



Mn jm 

dxD&f(xD)cx / Af^^j['{y)dy 



oc 



/ d(cos0)(l + i^^)P2(cos6 
J m 



= 0, 



(3.41) 



where the orthogonahty property of the Legendre polynomials is used. 

A deviation from zero, caused by the 0-functions is not large compared to 1, but anything is large compare 



to 0! One can artificially adjust formula (3.4C) to fulfill this sum rule. For instance, small corrections to the 
normalization of both terms with M = 1 and M = can be made to satisfy the sum rule in the form ( 3.32| ). 
However, the situation with the second sum rule, ( |3^ ) and (|3l^ , is more difficult and can not be fixed by 



any simple adjustments of the normalizations or ignoring the & functions. Similar to eq. (3.41), it can be 
written (neglecting the 6-function!) as: 



Md /m 



y cixD^f (xd) cx J yAf 



oc / d{cosO){MD ~ w 
-1 



(3.42) 



Thus, there is no reason for this sum rule to be satisfied with the nonrelativistic distribution function (3.40). 
Therefore the nonrelativistic formulae, in principle, violates the sum rules for the SFs hi2- However, we can 
still hope that it will be a small effect, one not noticeable in practice. 

Nonrelativistic formulae for other spin-dependent distribution, f^^^, are also a straightforward result 
of using density (2.16) instead of the relativistic density in the pole approximation in eq. ( [3.27 ) (see also 



footnote to the formula ( 2.19 )): 



f. 



'Tiy) 



5\y 



e(y) 



(3.43) 



w'(IpI) - l^'m + P2{cosd)w{\p\) (^w{\p\) - V2u{\p\) 



and the sum rule follows from eq. ( 2.1 
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Mo /m 



(3.44) 



F. Calculation of distributions: relativistic vs. nonrelativistic 



In order to understand if relativistic distribution functions, ( |3.21 ), ( 3.25 ) and ( ^j.27| ), are significantly 



different from the nonrelativistic distributions, (3.39), (3.4C) and (3.27), we have to understand fully an 



effect of the one pole approximation on the densities in the BS formalism. Indeed, discussion of Section [I B 
prepared us to the fact that we can not expect significant physical effects from the form of the densities, 
since both relativistic and nonrelativistic realistic models lead to similar results and we do not detect any 
special trend for the relativistic model. 

In Figs. ^ and ^ we compare the exact densities and densities in the one pole approximation, for two most 
important densities (70) and (7573). We find that one pole approximation leads to a significant change in 
the density. The densities in one pole approximation (solid curves) have much harder tail compared to the 
exact densities (dotted curves), starting at |p| ^ 0.5m, and effect is an order of magnitude at |p| ~ 1.5m. 
This can be qualitatively understood, considering an example of the charge density, (70). Indeed, picking the 
nucleon pole in the full integral, corresponds to neglecting the antinucleon (negative) contribution to total 
charge density which concentrates at high momenta, |p|. The presence of 0-functions in the expressions for 
the distribution functions, J^l^ and j^l^ , cuts off a part of the high momentum region, but this is a minor 
effect. Effect of "softening" caused by the 0-functions is also illustrated in Figs. || and ^ (dashed lines). 
These curves are the BS densities in the pole approximation (solid lines) multiplied by 0(Af£) — to +P3) 
integrated over cosQ. 

The results presented in Figs. |^ and pimply that the relativistic densities appearing in the formulae for the 
effective distribution functions are enhanced at medium and high momenta. Comparing these results with 
those in Figs. ^ and ^ we find that the effect is much larger than any model differences. Still, since the effect 
is concentrated at high momenta it is not clear if it leads to observable effects in deep inelastic scattering. 
To clarify this, in Fig. |l^ we present the effective distribution function, f^^^{y)- For completeness of the 
illustration, we compare distribution functions calculated using three different charge densities: (i) nonrela- 
tivistic density of the Bonn potential (dotted line), (ii) relativistic BS density in the one pole approximation 
(solid curve) and (iii) exact densities within the BS approach (dashed curve). The last curve is aimed to 
illustrate differences in the description of the mechanism of the reaction in nonrelativistic and relativistic 
approaches, i.e. it presents result for made up situation, if charge density of the nonrelativistic model is 
exactly the same as of the BS approach. We find that consistent relativistic description gives the effective 
distribution which is systematically harder at high momentum fraction, y. This is a result of harder charge 
density in the one pole approximation. It is interesting that the relativistic distribution is also enhanced at 
small ?/ — !■ 0. This effectively corresponds to larger "binding effects" in the BS approach which was observed 
in ref. j^. Very similar effects can be observed in other effective distribution functions, f^l^ and A/^/^. 
However, in these cases those effects are not so explicit, because of the oscillating nature of the distribution 
functions. (See e.g. discussion about A/^/^ below.) 
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IV. NUMERICAL CALCULATION OF STRUCTURE FUNCTIONS 



A. Unpolarized deuteron, 

The spin-independent SF of the deuteron, is calculated using the effective distribution functions 



presented in Fig. |1C| (see discussion at the end of the Section [II F ). The results of this calculation are shown 
in Fig. |l^ in the form of ratio of the SF of the deuteron and the nucleon. The nucleon SF, F^ is taken 
from ref. at — 10 GeV^. We find that BS approach gives a behavior of the deuteron SF qualitatively 
similar to the results of the nonrelativistic and "made up nonrelativistic" calculations. However, there are 
two delicate, but essential differences: 

• First, the ratio F^ / F^ in the BS approach is lower than others at small x, x < 0.5. This effect can 
be easily understood from the form of distribution, f^^^{y), in Fig. |l0| and formula ( 3.23 ). Indeed, at 
a; = the SF F^ is: 

Md /m 

F2^(0) = F2^(0). / r/^(y)dy. (4.1) 



Mo /m 

Therefore, F^{0)/F^{0) = 1 because / f^l'^{y)dy = 1 is a normalization, ( |3.30| ). When we move 





from X = to larger x, we "lose" part of the normalization, since the lower limit of integral in cq. (3.23) 
is X. Since the relativistic /^/^ is larger at small x than the nonrelativistic ones, it leads to a faster 
decrease of the relativistic SF with increasing x. The fact that the ratio F^ /F^ is less than 1 at small 
and intermediate x is know to be a result of the "binding of nucleons" (see e.g. p^-p7|,p| and references 
therein). 

• Second, the relativistic SF displays sharper rise at higher x, x > 0.5 than the nonrelativistic ones. 
This, again, can be understood from the form of the distribution. Fig. and the convolution formula. 



( 3.23 ). With increasing of x the role of the high momentum tail of the effective distribution is getting 



more important for integration in eq. (3.23) and at a; > 1.0 the tail is completely dominating. The 



deuteron SF at a; > 1 is presented in Fig. |l2|. 

It has been shown recently, that relativistic calculations lead to the larger binding effects than the non- 
relativistic ones P,^. However, the result of ref. joj which presents effects two to three times larger than 
ours still is neither explained nor physically understood. In ref. ||^ the size of the effect was not so large, 
but the method of numerical calculations was essentially based on the numerical inverse Mellin transform 
of the non-analytical function. The approximations had to be made, which eventually was not good at high 
momentum. That is why we specially attempted to verify quantitative size of the effect. It is worth to 
remember that we discuss quite tiny effect of ~ 1 — 2% in the ratio F^ /F^ . We also can measure the bind- 



ing effect using the energy-momentum sum rule (3.2E). The quantity Sn is a natural parameter controlling 



the binding in any particular calculations. For example, the nonrelativistic formulae allow for an analytical 
estimate of which essentially gives us an idea about size of the effect ]l6|,p7|p|: 
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where (T) is a nonrelativistic mean kinetic energy of nucleons in the deuteron. For the reahstic models 
typically (T) « 15 MeV, which gives (5 « 5 • 10^^. Calculating with the BS effective distribution of the 
present paper, we find S ~ 0.7 • 10^"^. Note, that in ref. ||^ the quantity (5 « 1 • 10^^ has been reported. 
We attribute this small discrepancy to the less sophisticated numerical approximation which had been made 
in 0]. This approximation slightly underestimates a high momentum behavior of f^/^ , 



B. Polarized deuteron, and 62 



The SFs of the deuteron (^) ^-^^ calculated within both the relativistic and nonrelativistic approaches. 
The relativistic calculations are based on the formulae ( 3.21 ), ( 3.22 )-( 3.25 ). The nonrelativistic calculations, 
eq. ( 3.40| ), uses the realistic wave function of the deuteron in the Bonn potential |l^. Another ingredient 
of the calculations, the nucleon SFs F(^2{^i Q^)i is again taken from ref. |Q at = 10 GeV^. The results 
are neither very sensitive to the particular choice of the parametrization for the nucleon SFs nor to the 
Q^-dependence of them. 

The distribution functions A/^/^(y) are calculated and the results are presented in Fig. |l^. The relativis- 
tic (solid line) and nonrelativistic (dotted line) calculations give similar behavior of the distribution function. 
Indeed, it is difficult to distinguish between them, not speaking about making definite conclusions. The third 
line in the Fig. |l^ is given for illustration, and presents yAf^^^{y) for the relativistic calculations. The 
calculation of the sum rules is more representative. To understand the scale of effects, which are discussed 
below, it is customary to define auxiliary quantities: 

Mo/m Mo/m 

J Abs (A/^/^(y)) ^ J Ahs (yAf''^^ (y)) dy ^ O.U. (4.3) 



The Bethe-Salpeter and nonrelativistic Bonn calculations give the same result in (4.3), with accuracy of 
^ 5%. Thus, the effective distribution functions, A/^/^, are an order of magnitude smaller than the usual 
spin-independent distributions f^/^ normalized to 1. This is not a very important circumstance, but it works 
against accuracy in the numerical calculations, since is a difference of two functions normalized to 1 

(M = 1 and M = 0). Numerically, the sum rule (3.34) (see also eq. (3.41)) is satisfied both in relativistic and 
nonrelativistic calculations with good accuracy, despite the approximate numerical "inverse Wick rotation" 
and the discussion after eq. ( ^.40 ). The corresponding integrals are ~ 5 • 10^'' and ^ 3 • 10^^ and they should 
be compared to the estimate (|4.3|). The sum rule ( ^.28 ) may be used to improve distributions A/^/^ by 
making integrals for ,fi^^ and f^^^ exactly the same. However, this does not lead to a significant variation 
of results for SFs, except x for fef (x). 

The behavior of 6f(a;) at a; ^ deserves to be considered more closely, especially for numerical calcu- 
lations, since the nucleon function {x) can be divergent at small x. Unfortunately it is impossible to 
estimate 6f (0) exactly for the realistic SF . However, a contribution of the singularity can be evaluated. 
Indeed, let us assume a singular behavior as Fi ^ C/x, then eq. (3.22) leads for small x to 
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Mo/r. 



Mo/r. 



C-Z 



-CA/^/-°(0) 



(4.4) 



where Z = in exact relativistic formula, but it can be a small number in numerical calculations or in the 
nonrelativistic formalism. Thus, the limit of the deuteron SF &f (x) as a; — *■ is a constant, but one has 
to exercise great care in performing numerical computations, since any error leads to a divergent behavior 



at small x. In this context, an adjustment of norms of the two terms in formulae (3.25) and (3.40) has a 
meaning of subtraction of the numerical error from &f at small x. 



The situation with the second sum rule (3.35) is quite different. Numerically it is violated more significantly 
than the previous one. Corresponding integrals arc ~ 1-10"'^ and ~ 3- 10^"^ for relativistic and nonrelativistic 



calculations respectively, i.e. about 0.7% and 2% compare to (O). Therefore, numerical approximations 
slightly damage the relativistic formula too. It is attributed to the numerical rotation in the Minkowski 
space. An adjustment of the normalization, as it has been discussed, slightly improves the accuracy (to 
0.5%). On the contrary, the result for the nonrelativistic approach is stable with respect to any adjustments, 
since it is defined by the formulae ( |3.42| ). 

The SFs bi and 6f are calculated within two approaches as well. The results are shown in Fig. |l^ a) 



and b). The behavior of the functions in Fig. ^ a) suggests the validity of the sum rule ( 3.34 ). At the 
same time, the nonrelativistic calculation for in Fig. |l^ b) (dotted line) obviously does not satisfy the 
sum rule dJI). The main difference of the relativistic and nonrelativistic calculations is at small x, where 
these approaches give different signs for the SFs. To illustrate the effect of the presence of the 0-function 



under integral in the nonrelativistic formula ( 3.40 ), the calculations have been done as well with a restricted 
interval of integration over |p|. The condition |p| < 0.7 GeV corresponds to the "softer" deuteron wave 



function, but makes the sum rule ( |3.42| ) exact. Corresponding SFs are shown in Fig. 14 a) and b) (dashed 
line). The result of this "experiment" is that the effect of 8-function is not quantitatively significant. It also 



does not affect the principle conclusion about the second sum rule (3.35), but makes the defect a tiny bit 



smaller. This is understandable, since the sum rule breaking terms in ( 3.42 ) is cx |p|cos0 



C. Polarized deuteron, 



The spin-dependent SF of the deuteron, , is calculated using the same three models as the spin- 
independent SF, F2: the fully relativistic BS approach (solid line), nonrelativistic approach based on the 
Bonn wave function(dottcd line) and nonrelativistic approach which uses the exact density of the BS approach 



(dashed line). The nucleon SF, gf' is taken from ref. |43|. The results of calculations are presented in Fig. |T^. 

For illustration we also present in Fig. |l^ the quantity (7573)f/^i which corresponds to a "model" for the 
deuteron SF (dot-dashed straight line ): 

g^x^Q^) = l.(^,^,)fl^.g^{x,Q'). (4.5) 
If we replace (7573)f/^i by the factor (1 — 3/2 • wd) from eq. (2.18), we get the formula usually used by 



experimentalists aiming to obtain the neutron SF, 5" from the combined proton and deuteron data. 
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Fig. lU shows that picture is different from the naive estimate ([i.5|). However, within the present day 
experimental errors it may be a reasonable approximation (see e.g. ref. |^). Huge leaps of the ratio around 
the constant (7573)^=1 at x < 0.7 are not too important. They correspond to zeros of the nucleon SF 
which are slightly shifted by the convolution formula. Systematic difference in the ratio exist between the 
nonrclativistic calculation (dotted line) and two calculations based on the BS densities (solid and dashed 
curves). This is a difference between D-wave admixture in the Bonn potential {wd ~ 4.3%) and our solution 
of the BS equation {wd ~ 5%). The rise of the ratio, at x higher than 0.7, is of the same nature as in the 
spin-independent case in Fig. nTl it is caused by the Fermi motion. 



V. SUMMARY 



In this paper, we have presented a study of the deep inelastic electron scattering on the deuteron in the 
Bethe-Salpeter formalism in the realistic meson-nucleon model. In particular, 

• The connection of the structure functions to the densities of the appropriate charges and currents is 
analyzed. By analysing the same densities in the nonrclativistic approach, we have systematically 
compared the relativistic and nonrclativistic calculations, and established sources of the relativistic 
effects. 

• Using our numerical solution of the Bethe-Salpeter equation amplitude with a realistic kernel, the 
leading twist structure functions of the deuteron, , and , are calculated in the fully relativistic 
fashion. 

• Our numerical calculations of the structure functions emphasize a qualitative agreement with previous 
non-relativistic results. However, we have found effects systematically distinguishing a consistent rela- 
tivistic approach from the nonrclativistic one: in the relativistic formalism (i) the magnitude of binding 
effects is larger and (ii) the effect of Fermi motion at high x is stronger; and (iii) the nonrclativistic 
calculations suffer unavoidable internal inconsistencies, which lead to small effects in the structure 
functions and gf , but seriously damage structure function and noticably break sum rules for 
this function. 

The present paper concludes our systematic study of the deep inelastic electron (muon) scattering on 
the deuteron in the Bethe-Salpeter formalism. The results are collected in the two papers, ref. and the 
present paper, and in part have been also published previously in refs. [p|-^|j3C|] . The main lesson, we have 
learned from this study, is that the deuteron in the deep inelastic reaction indeed behaves as a very slightly 
relativistic system. One has to look for special conditions or kinematics of the reaction to be able to find 
noticible relativistic effects. We have found certain situations where the relativistic approach is absolutely 
essential and use of the nonrclativistic methods is not justified. Most representative examples are the high 
X behaviour of the structure funtions and the spin-dependent structure fmictions, 6^*2 ■ 

Apart from the phcnomenological differences between the relativistic and nonrclativistic approaches, the 
most important merit of the covariant formalism is its consistency. Its close connection to the field theory 
guaranties that the calculated observables obey the sum rules and other general properties imposed by the 
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fundamental principles. In this sense the relativistic approach is definitely more advanced theoretically than 
its nonrelativistic counterpart and provides a better understanding of the deep inelastic scattering on the 
deuteron. 
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APPENDIX A: EXPLICIT FORMULAE FOR DENSITIES 

In this Appendix we present formulae which allow to restore explicit form of the various densities in the 
BS formalism. For convenience, we define auxiliary "densities" : 

27r 



} =^5l/^<^Tr{^'Mbo,p) I |*M(po,p)(P2-m)}, (Al) 



rN/D 




2tt 



, y"d0Tr{*M(po,p) I ^° j>*M(po,p)(p2-m)}^^^^-y"d0Tr{...}j,,,^o, (A2) 



I jJ/D^l I = /d</>Tr{§M(po,p) I Um(po,p)(P2 - m)},,_^ , (A3) 
[ fs (P) ) { [ 7573 J 

where integration over 4> leads to the trivial factor of 27r, since none of the matrix elements on the r.h.s. 
depends on 0. To obtain explicit expressions for the densities discussed in the paper we have to compare 



eq. (Al)-(A3) with definitions of corresponding densities. 

Note, that a formalism presented in refs. [Qj2| and in this paper can be easily adopted for an analytical 
computer calculations. The following results have been obtained utilizing the Mathematica package p^ : 

fo^^ip)^ ™ (-8V'ao(PO,P)V'to(PO,P) 8-002(^0, P)'0i2(PO,P)) (A4) 

I -8^piiPO:P)AoiPQ,p) /2 



\ ^/3 V3 

-V'ao(Po,p)V'i>i(Po,J3) H I 

{PO,P) +'2ll^ao{PO,P) +'2^a2{P0,p) +2lppi{pQ,p) 

+8V'to(po,p)^ + 8-0t2(po,p)^ + 8V'4°i(po,p)^ + 2ipvi{po,pY 
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/3^/^(p)= cos(^)| (A5) 
[^Md-po) \— 8y-V'pi(po,p)V't2(po,p) 

~ Vs Ts 1 

+ I — ^ +4y-V'a2(po,p)V'pi(Po,P) 

,0/2",/ ^ N , 8V't2(P0,P)V'i;l(P0,P)\ 

gf/'fobo.PjV't^ibcP) H I 

/ 0/ n2 2V'a0(P0,P)^ 8\/2Va0(P0,p)V'a2(p0,P) 
+ P \^i^al(PO,P) ^ ^ 

, 2V'o2(P0,P)^ „ , . ^2 , 8'!/'io(PO,P)^ „ , . ^2 

+ ^ ^ -ZV-pilPcP) + 2ip^i{po,p) 

32s/2tpto{po,p)ipt2ipo,p) Sipt2iPo,pf ,0,0/ \2 

+ o o — + ^y^ti {po , p) 



/^/^(p) = com I (A6) 



+ m 

+ 



(^-4y6f/'to(po,p)'0t>i(Po,p) - 4V3-(/'t2(po,p)'i/'«i(Po,p) 

P (2 V'aO (PO , P) ^ + 2 \/2 VaO (PO , V'a2 (PO , P) + ^"02 (PO , P) ^ 

-8tpto{Po,pf - 8\/2Vto(Po,p)V't2(Po,p) - ^■4>t2{po,pf 
-12V't°i(Po,p)^ + 3V'«i(po,p)^ 



fT''ip)= P2(cos(0))i (A7) 



- 4:\ -lljpl(j>0,p)'lpt2{P0,P) 



+ m 



\ Vs V 3 

r, , ( ^1 I X 8V'a2(P0,P)V'^l(P0,P) 
-2Y 3V'a0(P0,p)V'i;l(P0,p) 

(^4\/2 Va2 (po , P) V'tO (pO , P) + 4\/2 VaO (pO , p) V't 2 (PO , p) 
+8^i'a2(P0,p)V't2(P0,p) + 4\/2V'°i(po,p)V'?i(PO,P) 

-2\/2il}px{po,p)il)yi (po, p)) 

+ (|Md -po) (-2V^Va0(P0,p)V'a2(p0,P) - 2lpa2{po,pf 
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+ p 



-8\/2ipto{po,p)'^t2{po,p) - 8ijjt2{po,pf 

-&1ptl{P0,pf -^tpvl{PO,pf^ I 



l8'^pi{po,p)i^to{po,p) , . /2 , / ^ ^ 



A ^ I I Ml ^ , 2V'a2(P0,P)V'j;l(P0,P) 



A/o"/"(p) = 



+ m {8tpao{PO,P)tlJto{PO,P) - ^■4'a2{P0,P)'^t2{P0,P) 

-4\/2V'°i(po,p)V'n(Po,p) + 2V2tppi{po,p)MPo,p)) 

+ i^Md-po) (^-2tpao{po,pf +1pa2iP0,pf -StptoiPO,pf 

+^'>Pt2{po,pf - 4V'°i(po,p)^ - '>PviiPo,pf^ 

P2(C0S(^)) I 

(^-12V2yja2iPo,p)^to{po,p) - ao{PO,P)'ipt2{po,p) 

+ l2lpa2{po,P)lpt2{PO,P)j 

(8V3tppi{po,p)tpto{po,p) - &V6ilJpiipo,p)->Pt2{po,p) 

+2\/6tpa0{P0,P)'>Pvl{P0,P) + 2\/3Va2(P0,P)V'«l(P0,P)) 
ao{PO,P)i'a2{PO,p) " 3^)02 (PO,P) 

-eV-pife),!?)^ + 24\/2V'tobo,p)^t2(po,p) - ^2tpt2{Po,py 

n 2 2^ 



m 



+ P 



+ 



cos(^) 

P (^-4'^a0(P0,p)^ + 2V2lpa0ip0,P)lpa2iP0,p) + ^tpa2{po,pf 

+16'ipto{Po,pf - SV2ipto{po,p)tpt2iPo,p) - 16tpt2{po,pf^ 
+ i^Md-po) (-8V3Vpi(Po,p)'0to(po,p) - ^VQil^pi{j)o,p)'4>t2{po,p) 

-2V6lpa0{P0,P)'>Pvl{P0,P) + 4:VSlpa2{P0,P)lpvl{P0,P)^ 

^4^3^00 (PO,P)V'pl(PO,P) + 2\/6ilJa2iPo,p)ippi{po,p) 
+'iV6iptoiPo,p)'4'viiPo,p) - S>V^1pt2{P0,P)''Pvl{P0,P)j 



+ m 



+P2(COS(0)) 



22 



+ m (^-6V6->pa2iPo,p)tppi{po,p) + l'2V3ipt2{po,p)'>Pvi{po,p) 

+ P (^-^^aliPO^pf - 9-0a2(PO,P)^ + G^pliPO,pf 

+S6'>pt2{po,pf + 12-0t"i(Po,p)^ - ^■ipviiPo.pf ) 



APPENDIX B: HADRONIC TENSOR FOR THE DEUTERON. PROJECTORS. 



The parametrization of hadron tensor for the deuteron, utihzed in the present paper, is given by eq. (3.1). 
It has both symmetric, {..•}, and antisymmetric [. . .] parts in respect to the permutation of its indices: 

Three physical vectors are used in this parametrization: 

1. Pjj the deuteron momentum. In the rest frame of the deuteron Pjj — {Mjj, 0); 

2. q the momentum transfer in the deep inelastic scattering. With proper choice of the orientation 
of the coordinate system q = (i^, 0, 0, — + Q^)- In the deep inelastic limit, when jv^ 0, 
pq = j/(po +P3)- 

3. SniM) is the total angular momentum of the deuteron, i.e. spin of the deuteron as an elementary 
particle: 

SUM) = -j^e'^f'^'E;iM)E^{M)PD5, (B2) 

f ^(0,-l,-*,0), Af = l 
E{M) = } (0,0,0,1), Af = (B3) 

[ ^(0,1, -z,0), M = -l 

The symmetric part of the hadron tensor, W^^^^ , contains terms proportional to the two following tensor 
structures: 

Tl.])--9,.+ '-f'. (B4) 

Because of the conservation of the electromagnetic current, the contraction of the hadron tensor with in 
respect of any index is zero. That is why only PdiiPdh and g^^i, are available to construct the projection 
operators to extract structure functions F-^^ and from the hadron tensor. We introduce the following 
coefficients: 
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Dc = Ci -04-02- O3, (B7) 



'21. 

d2 

Then the structure function are recovered by: 



A^^g^'^W^^, A^^^^W^^. (B8) 



FF = (B9) 
FF = (BIO) 
The antisymmetric part of the hadron tensor of the deuteron, W^^^ , in the general case has a form: 

^[H = ^^-"Z'^" {^oiM) {gn^D,Q')+g?{xn,Q')) - P'n^-^^^9?i^D,Q')] , (Bll) 

where the second structure function, g2, is vanishing in the deep inelastic limit, v ^ 00, — > 
00, Q^/i^ — > const. We do not discuss this structure function in the present paper. To obtain the spin- 
dependent structure function, , we construct the antisymmetric projectors. The two following projectors 
are equivalent for our purpose: 

= te.^^pq'^S^^iM), (B12) 



^ i{SD{M)q) 

It is interesting that in the limit jv^ — > 



%u = (^ixvaffQ t^D- K^^-^) 



D ^ Mfl ^ 
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meson 


coupling constcints 


mass 


cut-off 


isospin 
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OD/f47r): \q„IqA 

^ B / V / ' \_>JV / nil 


ur. GcV 


A.GcV 






12.2 


0.571 


1.29 





S 


1.6 


0.961 


1.29 


1 


TT 


14.5 


0.139 


1.29 


1 


V 


4.5 


0.549 


1.29 







27.0; [0] 


0.783 


1.29 





P 


1.0: [()] 


0.764 


1.29 


1 




;// = 0.9.39 GcV, zd = -2.21 


5 A1(A' 





Table 1. The set of parameters of the kernel of BS equation used in this work. 
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FIG. 1. Triangle graph for nucleon contribution to the deuteron matrix element of operator O. 



FIG. 2. The charge density in the deuteron, eqs. (2.9) and (2.11) calculated in different models, nonrelativistic and 
relativistic. Curves: the Bethe-Salpeter amplitude (solid), the Bonn wave function (dashed), the Paris wave function 
(dotted). 

FIG. 3. The (73)(p) in the deuteron calculated in different models. Curves; BS amplitude (solid), the Bonn wave 
function (dashed), the Paris wave function (dotted). 

FIG. 4. The "tensor density" of the deuteron, {7o)m=i(p) - (7o)m=o 

(p), calculated in different models. To exclude 
the angular dependence the densities are divided by P2{cos6). Curves: BS amplitude (solid), the Bonn wave function 
(dashed), the Paris wave function (dotted). 



FIG. 5. The spin density, {7573)a/=i(p), in the deuteron calculated in different models (see Section [IE for defi- 
nitions). To exclude the angular dependence for the present figure, densities are integrated over all angles. Curves: 
BS amplitude (solid), the Bonn wave function (dashed), the Paris wave function (dotted). 

FIG. 6. The one-pole contribution to the charge density of the deuteron, (7o)^ofe(p)i calculated with the BS 
amplitude. Curves: the leading term in the Taylor expansion (dotted), first two terms, up to oc p2 (dashed), first 
three terms, up to oc pi (solid) . 

FIG. 7. The one-pole contribution to the spin density, of the deuteron, {7573)^ofe(p)5 calculated with the BS 
amplitude. Curves: the leading term in the Taylor expansion (dotted), first two terms, up to oc (dashed), first 
three terms, up to oc pi (solid). 

FIG. 8. Different "versions" of the charge density of the deuteron calculated with the BS amplitude. Curves: 
exact, {7o)^'^(p) (dotted); in the pole approximations, (7o)^ofe(p)5 (solid). Dashed curve presents effective density 
for deep inelastic scattering which is the density in the pole approximation with account for the O-function cut-off 
at high momenta. 



FIG. 9. The same as in Fig. |8| but for the spin-density, (7573)^'^(p). 

FIG. 10. The effective distribution for the nucleon contribution in the deuteron structure function , /^''^(y). 
Curves: fully relativistic BS (solid); nonrelativistic Bonn (dotted), nonrelativistic, but using the BS charge density 
(dashed) . 
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FIG. 11. The ratio of the deuteron and nucleon structure functions, / F2 , calculated in different models. Curves 
correspond to three effective distributions from Fig. |l^: fully relativistic BS (solid); nonrelativistic Bonn (dotted), 
nonrelativistic, but using the BS charge density (dashed) . 

FIG. 12. The deuteron structure function {x) at large x, calculated in different models. Curves correspond to 
three effective distributions from Fig. |lo| fully relativistic BS (solid); nonrelativistic Bonn (dotted), nonrelativistic, 
but using the BS charge density (dashed). Dash-dotted curve present the free nucleon structure function, F2 ■ 

FIG. 13. The effective distribution functions for the deuteron structure functions, ^/^^°(j/)i calculated in 
different models. Curves: fully relativistic BS (solid); nonrelativistic Bonn (dotted). The dashed curve presents 
J/A/^/^(j/). 

FIG. 14. The deuteron structure functions bf (a) and 62" (b), calculated in different models. Curves: fully 
relativistic BS (solid); nonrelativistic Bonn (dashed). The dotted curve presents 6^2 calculated with the "soft" 
nucleon distribution, the Bonn distribution but cut-off at |p| > 0.7 GeV . 

FIG. 15. The ratio of the deuteron and nucleon spin-dependent structure functions, g\ /g\ , calculated in different 
models. Curves correspond to three effective distributions: fully relativistic BS (solid); nonrelativistic Bonn (dotted), 
nonrelativistic, but using the BS spin density (dashed). 
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